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By 
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Io  General  Discussion 


Between  engineering  analysis  and  flight  testing  there  are  several  tests 
which  help  to  establish  the  flutter  characteristics  of  the  given  configuration* 

The  tests  may  be  of  a  static  or  dynamic  nature,  and  some  systems  are  so  complex 
that  none  of  the  tests  define  the  structure  well  enough  to  show  its  flutter  tend¬ 
encies,  It  is  not  the  purpose  of  this  report  to  describe  the  various  possibilities, 
but  to  propose  an  additional  technique  for  ground  testing  by  a  simulation  of  flutter* 

It  has  been  learned  that  in  certain  cases  the  distributed  aerodynamic 
forces  on  a  surface  may  be  represented  by  a  single  concentrated  force  acting  180 
degrees  out  of  phase  with  the  angle  of  attack*  The  flutter  of  a  two-dimensional 
supersonic  tail  heavy  surface  can  be  understood  on  this  basis,  and  it  is  implied 
that  the  structurally  three-dimensional  case  can  be  handled  equally  well*  A 
vibration  test  of  the  actual  surface  can  be  used  to  establish  the  magnitude  of 
the  "simulator  speed" „  It  is  one  step  further  to  use  a  feed-back  system  to  pro¬ 
duce  a  vibration  representative  of  the  same  flutter  mode.  It  involves  using  an 
electrical  signal  representing  the  angle  of  attack  to  produce  a  shaker  force  180 
degrees  out  of  phase  with  the  signal *  When  the  gain  is  raised  to  a  certain  critical 
level,  a  self -excitation  is  produced.  The  corresponding  shaker  force  and  the  sur¬ 
face  displacements  can  then  be  used  for  a  direct  computation  of  the  flutter  speed. 

The  technique  is  currently  being  applied  to  a  simple  flat  plate  cantilever  for 
which  the  flutter  speed  can  be  computed;  the  comparison  of  the  simulator  speed 
and  the  computed  speed  is  needed  to  show  that  the  device  is  operating  satisfactorily* 

Nore  adequate  simulation  can  be  achieved  by  the  use  of  more  pickups  and 
more  shakers,  so  that  the  forces  and  moments  produced  will  correspond  more  accu¬ 
rately  to  those  produced  by  an  oscillating  airfoil* 

II*  Analysis  Using  Computed  Data  for  Uncoupled  Modes 

It  is  wefll  to  begin  with  a  simple  three-dimensional  surface  having 
bending  and  rotational  degrees  of  freedom,  reacting  to  an  aerodynamic  force  at 
the  elastic  axis.  Figure  1. 


FIGURE  I 


units  msics  uimrwr 
m  itm  nonius  wsnr'«iit 
sunn  imu  ktiruio 


2. 


The  application  of  Lagrange* s  Equations  produces  the  tw:  equations 
required  for  equilibrium: 


(1)  (kfr  -  Mq  co2)  q_^  -  S  o>3q2  «  /  p  Di  c  dy 

■  ~N  Qq 

(2)  -  Sco2qi  +  (ko  -  I®2)  q2  ■  0 


For  a  distributed  spring  mass  system  these  quantitie )  have  the  nature 
of  generalised  forces,  defined  on  page  io# 

A  condition  of  simulated  flutter  exists  when  the  de  ..errainant  is  equal 

to  zero. 

(3)  kh  kj^-  Gh  Ga  -  <o2^M  kao+  I  k^j  +  M  I  a1 

+  S  o^L  -  S2  a>4  +  i  Hco*]+  i  Gh  I  co2]  =0 

The  real  and  imaginary  parts  are  separately  set  equal  to  zero: 
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This  value  of  the  frequency  is  then  used  in  the  real  part  of  the 

equation: 
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The  equation  ultimately  simplifies  to: 
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The  equation  can  be  put  into  a  form  suitable  for  correlating  results 
with  NACA  calculations  by  using  the  definitions  on  page  *£?-  and  the  ultimate 
result  is: 
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3.51 
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7.85U 

3.90 
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*  uncoupled  bending  frequency 

°  uncoupled  torsion  frequency 

b  «*  surface  semi-chord 

c  ■  chord 

x  «*  distance  of  center  of  gravity  aft  of  the  elastic  axis 
r  «*  radius  of  gyration  about  elastic  axis 

M-  n  m/lj  b^  p 

m  *  mass  per  unit  span 

P  a  density  of  air 

M  "  Mach  number 

n 

Note:  The  elastic  axis  is  at  mid-chord  in  this  example. 
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III,  Analysis  Based  Upon  Measured  Data  from  Forced  Vibration 

It  has  been  shown  that  the  two  dimensional  simulator  speed  has  a 
correlation  with  the  NACA  data.  The  next  step  is  to  show  that  experimental 
data  can  be  obtained  from  a  distributed  spring-mass  system  which  can  show  that 
the  real  system  is  in  a  condition  of  simulated  flutter.  This  involves  the 
application  of  a  single  shaker  force  at  some  point  on  the  elastic  axis  of  the 
surface,  with  the  assumption  that  the  elastic  axis  and  the  centers  of  pressure 
coincide.  The  angular  response  of  the  system  can  then  be  monitored  to  find 
the  condition  wherein  the  angular  response  is  180  degrees  out  of  phase  with 
the  applied  force. 

If  the  system  vibrates  in  this  particular  shape  at  some  supersonic 
Koch  number,  the  distributed  aerodynamic  forces  produce  a  generalized  aero¬ 
dynamic  force.  And  if  the  force  thus  developed  is  equal  to  the  generalized 
force  produced  by  the  single  shaker,  it  is  considered  that  such  forces  may 
effectively  replace  the  shaker  as  a  source  of  the  observed  vibration.  In  order 
to  verify  this  pointy  one  can  compute  these  forces,  equate  them.,  and  establish 
that  their  equality  is  true  if  the  original  determinant  is  equal  to  zero0  The 
following  equations  show  this  result. 

The  response  of  a  two  degree  of  freedom  system  to  a  force  at  the 
elastic  axis  is  determined  from  LaGrange's  Equations, 


(10)  (kh  -  Mq  to2)  q^  -  S  co2  q^  «  Pa  Dia 

(11)  -  S  co2  a  +  (k  -  I  co2 )  q  **  0 

j.  Ci  2 


The  response  of  the  bending  coordinate,  qit  is : 
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The  displacement  at  the  elastic  axis  is: 
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The  angle  at  any  spanwise  station  is; 
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A  condition  of  flutter  exists  when  these  two  forces  are  equal, 
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This  expression  is  to  be  compared  with  the 
by  an  expansion  of  the  flutter  determinant: 
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IV,  Flutter  Simulation 
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in  a  simulator  flutter  mode  exactly  the  same  as  the  one  reached  by  the  forced 
vibration  method.  In  this  case,  the  end  result  is  reached  more  quickly,,  A 
3t.udy  of  such  a  system  is  a  necessary  part  in  the  development  of  more  accurate 
flutter  simulation  by  the  use  of  many  shakers  and  pickups. 

The  amplifier  necessary  for  the  conversion  of  signal  voltage  to  shaker 
current  is  not  commercially  available.  Usually  a  phase  shift  is  present  between 
input  and  output  which  is  a  function  both  of  load  and  frequency.  Such  a  phase 
shift  has  been  eliminated  in  a  transistorized  amplifier  designed  by  George  Bush 
of  BTR/APL.  The  amplifier  is  approximately  flat  from  D.C,  to  1  KC,  and  prelimi¬ 
nary  tests  have  shown  satisfactory  phase  shift  characteristics.  The  development 
of  the  amplifier  is  the  key  to  the  whole  problem  of  the  simple  simulator.  For 
details  on  the  circuit,  reference  (a)*  should  be  consulted.  The  amplifier  input 
voltage  proportional  to  surface  angular  motion  may  be  produced  by  a  device  which 
measures  the  fore  and  aft  motion  of  the  tip  of  an  arm  mounted  perpendicular  to 
the  surface.  Method  A,  Figure  2,  This  is  one  simple  way  of  measuring  the  angle, 
but  it  is  valid  only  when  the  plate  is  not  deflecting  fore  and  aft0  A  recent 
inprovement,  Method  B,  Figure  (2),  has  been  to  measure  the  differential  output 
from  two  pickups  located  at  the  leading  •'nd  trailing  edges.  The  pickup  system 
is  comprised  of  a  Schaevitz  linear  differential  transformer  and  a  phase  comparator. 
The  comparator  output  voltage  is  used  to  drive  the  Bush  Amplifier.  The  amplifier 
output  voltage  is  used  to  drive  a  Goodjnan  V-U7  shaker  as  shown  in  Figure  3. 

The  operation  of  the  simulator  involves  a  slow  increase  in  the  gain 
of  the  system  until  the  aerodynamic  surface  oscillates.  The  force  acting  may  be 
found  from  the  appropriate  calibration  factor  and  the  shaker  current,  A  measure¬ 
ment  of  the  displacement  at  the  shaking  point  then  enables  one  to  compute  the 
generalized  shaker  input  force. 


(2^)  0_  =  o  h 

o  S  a 

The  next  step  is  to  compute  the  generalized  aerodynamic  force  that 
would  be .produced  by  the  surface  undergoing  the  same  displacements  at  some 
chosen  M^ch  number. 


L 

(2H)  Qa  °  Ct  q1  /  h  a  c  dy 

a  o 

The  expression  for  takes  on  an  interesting  form  when  the  centers 
of  pressure  lie  along  the  mid-chord,  and  the  leading  and  trailing  edge  dis¬ 
placements  hf  and  hj.  have  been  measured; 

(25)  «A  -  C^q1 


If  the  generalized  shaker  force  is  less  than  the  generalized  aerody¬ 
namic  force  a  condition  of  flutter  can  be  expected  at  the  chosen  Mach  number. 


^Reference  (a)  APL/BBE  Dwg.  57020 
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FIGURE  2 


Ground  Flutter  Simulator 
Schematic  Diagram 


Method  (A)  for  Sensing  Angular  Motion 
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As  the  simple  simulator  now  stands  it  is  a  tool  for  investigation  of 
various  simple  surfaces  to  verify  the  ideas  presented  here0  In  order  to  detect 
flutter  speeds  of  missile  surfaces,  the  sensitivity  has  to  be  greatly  increased 
by  the  use  cf  an  appropriate  voltage  amplifier 0  It  should  be  pointed  out  that 
the  simulator  has  a  wider  application  to  autopilot  stability  studies,,  One 
envisions  two  simulators  attached  to  symmetrical  wings  activated  with  flight 
servo,  gyro  and  allied  electronic s.  The  aeroelastic  high  frequency  instabili¬ 
ties  can  thus  be  investigated 0 

It  has  been  pointed  out  that  the  simple  simulator  is  a  tool  for  in¬ 
vestigations,  and  that  the  results  need  careful  evaluation  in  every  casec  In 
order  to  provide  more  adequate  force  and  moment  simulation,  more  shakers  and 
pickups,  together  with  appropriate  circuitry  for  producing  phase  shifts s  need 
to  be  added  to  the  system,, 

It  is  believed  that  use  of  the  simulator  will  be  helpful  in  the 
determination  of  the  flutter  speeds  of  missile  surfaces,  and  in  evaluating  the 
aeroelastic  stability  characteristics  of  assembled  missiles. 
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SYMBOLS 


10. 


kjj  «»  generalized  complex  stiffness  appropriate  to  the  uncoupled  bending  mode. 
kh  -  kh0  (1  +  i  gh) 

kh0  "  ^  <4 

gh  ■  damping  parameter 

kg  »  same  as  above  except  that  it  applies  to  the  uncoupled  torsion  mode. 

L 

Hr  “  C  /  ip  \  dy  3  <4 

o  o  v 

Mq  «  generalized  inertia  force  for  the  uncoupled  bending  mode. 

M  -  /L  m  dy 
8  o 

S  »  generalized  inertia  coupling  (measure  of  unbalance)  between  the  two  modes. 

S  «  /^ixD  d  dy 

I  *  generalized  inertia  element  for  the  uncoupled  torsion  mod8. 


*  amplitude  of  time  dependent  function  fo’-  the  bending  mode. 
q2  «  same  for  the  torsion  mode. 

P&  ■  amplitude  of  oscillating  force  at  frequency  applied  at  point  "a”. 

D  a  “  modal  displacement  at  point  "a"  due  to  bending,  (Note  that  "a”  is  on  the 
axis  of  twist.) 

Dx  **  bending  modal  displacement  at  any  point. 

D_  **  uncoupled  torsional  modal  displacement, 

u 

^La  "  supersonic  lift  coefficient,  yj  ~ 2  '£ 

q1  “  dynamic  pressure,  ^pv2 
P  ■  air  density, 

v  a  air  speed, 

C  a  total  chord  at  a  given  station 

y  a  spanwlse  coordinate, 

p  *»  pressure  acting  on  the  aerodynamic  surface 


